A simple model for some unusual properties of martensitic transformation 
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We report a detailed numerical investigation of a recently introduced two dimensional model for 
square-to-rectangle martensitic transformation that explains several unusual features of the marten- 
sitic transformation. This model includes inertial effects, dissipation, long-range interaction between 
the transformed domains and an inhomogeneous stress field to describe the effect of lattice defects 
which serves as nucleation centers. Both single-site nucleation and multi-site nucleation has been 
studied for single quench situation and thermal cycling. The final stage morphologies of single-site 
nucleation and multi-site nucleation bear considerable similarity suggesting that the initial dis- 
tribution of the defects is not important. Thermal cycling using continuous cooling and heating 
simulations show the existence of hysteresis in the transformation. More importantly, the rate of 
energy dissipated occurs in the forms of bursts with power law statistics for their amplitudes and 
durations which explains the results of acoustic emission signals observed in experiments. When the 
system is cycled repeatedly in a restricted domain of temperatures, the dissipated bursts of energy 
are repetitive, a feature observed in experiments. The associated morphology shows a complete 
reversal of the martensite domains thus throwing light on the mechanism underlying the shape 
memory effect. The model also exhibits tweed like patterns. 
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PACS numbers: 81.30.Kf,05.40.-a,64.60.Ht,45.70.Ht 



I. INTRODUCTION 



Motivation 



Martensitic transformation often exhibits unusual fea- 
tures that are not expected of a first order transition. 
One such effect is the well documented pretransitional 
effect observed as the system approaches the martensitic 
transformation temperature pi 0, 0, Q . The associated 
enhanced levels of fluctuation observed in several mea- 
surable quantities ( such as the anomalous scattering and 
tweed structure) has been recorded in a number of dif- 
ferent systems 0, 0, Q- This feature, however, is the 
signature of critical fluctuations in a second order transi- 
tion. As another example, consider an observation that 
relates to acoustic emission (AE), a feature that usu- 
ally accompanies martensite transformation. Recently, 
in experiments on Cu-Al-Zn single crystals, Vives et al 
reported that the distributions of the amplitudes of the 
AE signals and their durations obey power law statistics 
when the samples were subjected to slow thermal cycling 
[H, both during cooling and heating runs. The statis- 
tics were reported to be robust over a range of cooling 
and heating rates. Again, power law statistics imply scale 
free nature of the underlying process. However, in this 
case, as the AE signals are accumulated during thermal 
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cycling, there is no tuning as in second order transition. 
Thus, it is actually reminiscent of self-organized critical- 
ity introduced by Bak et al §. Another unusual and 
interesting property of the AE signals reported is the 
high degree of reproducibility and statistical correlation 
in time when the system is subjected to repeated thermal 
cycling over a restricted range of temperatures 9] . The 
near repetitive nature of AE signals during successive cy- 
cles have been shown to be correlated with the growth 
and shrinkage of martensite domains. Thus, exploring 
this correlated behavior should help us to understand 
the shape memory effect as well. Moreover, at a concep- 
tual level, the nature of correlation of the repetitive AE 
signals is significantly different from the power law na- 
ture of the statistics observed during full thermal cycling 
mentioned above. Thus, it would be desirable to capture 
the seemingly conflicting properties using a single model. 
The purpose of the paper is to report a detailed study 
of a recently introduced two dimensional phenomenolog- 
ical model that is shown to explain both the power law 
statistics and correlated behavior of AE signals [Tol . 
As we shall show, the model also captures the precursor 
effect. 



B. Background 

The martensitic transformations are technologically 
important class of first order, solid-solid, diffusionless 
structural phase transformations. Unlike other phase 
transformations where diffusion disperses the neighbour- 
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ing atoms, here, neighbouring atoms in the parent phase 
remain so in the product phase also [12|. However, the 
lattice gets distorted due to spontaneous displacement of 
the atoms ( from their positions in the parent lattice) ac- 
companying the discontinuous change in the shape and 
symmetry of the unit cell. This creates long-range strain 
fields which in turn strongly depend on the relative po- 
sitions and orientations of the martensitic plates. (Typ- 
ically, the martensite morphology consists of thin plate- 
like domains with twinned structure that are oriented 
along the elastically favorable habit plane directions.) 
Thus, the transformation path depends on continuously 
evolving configuration dependent long- range strain fields 
that eventually leaves the system in a metastable state. 
Thus, given a fixed quench, and hence a fixed amount 
of drive, the amount of transformed phase is fixed and 
further undercooling would be required to increase the 
product phase. For the same reason, the transformation 
occurs over a wide range of temperatures. On cooling, 
the transformation starts at a temperature M s , called 
the martensite start temperature and is completed at 
a temperature Mf (martensite finish temperature). In 
the reverse heating cycle, the transformation is initiated 
at a temperature A s (austenite start temperature ) and 
ends at a temperature Af (austenite finish temperature), 
which in general can be much higher than the martensitic 
start temperature M s . The related shape memory ef- 
fect usually accompanies the transformation. The above 
features also imply that thermal fluctuations have lit- 
tle role in the transformation kinetics. Such martensites 
are called athermal. It must be mentioned that marten- 
sitic transformations can be induced by applying external 
stress. The combined use of temperature and stress have 
found many practical applications |l3j| . 

In athermal martensites, the martensitic phase nucle- 
ates from isolated regions in the crystal which are usually 
defects like dislocations or grain boundaries 0, 0, ITq . 
Thus, quenched disorder plays an important role in the 
initial kinetics of athermal martensites. Another interest- 
ing feature of these transformations is the autocatalytic 
cooperative nature of nucleation 14] . It is observed that 
a martensite domain that has nucleated from a specific 
site in the crystal, triggers the nucleation of other do- 
mains in the vicinity. Such a correlated nucleation was 
observed by Ferraglio and Mukherjee who studied 
heterogeneous nucleation in Au-Cd alloys. 

A large bod y of theoretical work has accumulated over 
the years 0, TTij lvf\ . ( See for a recent summary 
of the status of theory.) Recent approaches use a contin- 
uum field theoretic models with strain or displacement 
fields as order parameters 0, OS, 0, • Within this 
framework, Cao et al \2\\ have dealt with the problem 
of heterogeneous nucleation at localized defect sites by 
simulating the influence of defects through an inhomogc- 
neous stress field. Recently, Wang and Khachaturyan |23 | 
have investigated the dynamics of improper martensitic 
transformations based on a time dependent Ginzburg- 
Landau (TGDL) approach. This approach has also been 



extended to proper martenstic transformation |24|. 

Here, we recall a few models relevant for the present 
work. The tweed structure has been explained as arising 
due to inhomogeneous distribution of the components of 
the martensitic alloy, for instance, Pd in Fe [22| using a 
strain based model. On the other hand, the power law be- 
haviour of avalanches occurring in first order transitions 
has been modelled by using disorder based Ising models 
[25j . In this case, the power law statistics of avalanches 
arises in the presence of a critical amount of disorder. 
Since quenched-in-disorder (defects) plays an essential 
role in the nucleation process of martensitic transforma- 
tion, it appears that these kinds of models [25| may be 
relevant to martensitic transformation. However, by sub- 
jecting the system to repeated hysteresis cycles, Vives et 
al, have verified that the system evolves towards a criti- 
cal state independent of the initial treatment of the alloy 
suggesting a dynamical evolution of the system towards 
such a state. According to these authors ||, the inter- 
pretation is that in real martensites, although there are 
quenched- in defects, the disorder that is responsible for 
the power law is actually generated during the transfor- 
mation itself. Indeed, this feature with lack of tuning 
is reminiscent of self-organized criticality, a word coined 
for slowly driven spatially extended systems evolving to 
a critical state where the statistics of avalanches obey a 
power law [8(. Since the introduction of this concept, 
there are large number of reports of physical systems ex- 
hibiting SOC like features, for example, earthquakes (2(|, 
acoustic emission from volcanic rocks |27I. stress drops 
during the Portevin Le-Chatelier effect |28| and biologi- 
cal evolution |2j| , to name a few. 

Recently, a simple two dimensional phenomcnological 
model has been shown to capture the power law statis- 
tics under thermal cycling This model attempts to 
incorporate the essential features of systems evolving to 
SOC state, namely, slow driving, threshold dynamics, ap- 
propriate relaxational mode without any recourse to tun- 
ing any relevant parameter. Surprisingly, this model also 
captures the correlated nature of the AE signals when 
the system is cycled in a restricted range of tempera- 
tures with the associated growth and decay of the 
martensite domains. 

The purpose of this paper is to present the results of an 
extensive numerical simulations on this two dimensional 
model describing square-to-rectangle martensitic trans- 
formation |l(J, [ll|. We discuss the results of a single- 
defect quench, multi-defect quench, thermal cycling over 
broad range of temperatures, the power law statistics 
arising in both in single and multi-defect cases, the cor- 
related nature of AE signals and the associated shape 
memory effect, and finally the tweed like structure. 

The organization of this paper is as follows. In sec- 
tion II, we describe our model starting from a free-energy 
functional and derive an equation of motion for a strain 
order parameter. The model has three free parameter 
one of which is the temperature. Section III contains the 
results of a detailed numerical simulations carried out for 
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various conditions of quench parameters. In section IV, 
we consider thermal cycling for single and multi-site cases 
both of which exhibit thermal hystersis. By calculating 
the energy dissipated during the transformation, we show 
that the distribution of the amplitudes and durations of 
the energy bursts during thermal cycling obey power law 
statistics as in experiments [6j. Section V deals with the 
correlated behavior of AE signals and its correspondence 
with shape memory features. Section VI, deals with the 
tweed structure. We end the paper with some observa- 
tions on the model. 



II. 



THE MODEL 



The basic idea of the model is to include all the im- 
portant features of athermal martentites such as inertial 
effects, long- range interaction and dissipation. The in- 
ertial effect is included by accounting for finite propaga- 
tion time in a manner similar to that considered by Bales 
and Gooding [3(], H3 • These authors have demonstrated 
that inertial effects prohibit the growth of martensite as 
a single variant in the presence of dissipation (in one 
dimension). Instead, the martensite grows as an alter- 
nating arrangement of the two variants. The reason for 
including dissipation stems from the recognition that the 
parent-product interface moves at near velocity of sound 
as suggested by the emission of AE signals. Associated 
with this movement, there is a dissipation which tends 
to relax the system towards local equilibrium. We in- 
clude this through a Rayleigh dissipation functional (3^| . 
( Note that conventional phase transformation take place 
at sufficiently slow pace thereby providing adequate time 
for quasi-steady state conditions to be attained in a short 
time.) We also include heterogenous nucleation at defect 
sites by including an appropriate strain energy [2l| . 

In order to capture these seemingly different types of 
features within the scope of a single model, it is desirable 
to keep the model as simple as possible. To this end we 
first consider a 2d square-to-rectangle transition [T^. l2fj| . 
However, the free energy in general depends on all the 
three strain components defined by 



with 



ei 



(?7ii + 7722) 

(Vn - 7722) 

V2 
V12 = 7721, 



dxi 



) 



(1) 



(2) 



referring to the components of the linearized strain ten- 
sor and tij's are the displacement fields in the direction 
i ( i — x,y). The components e\, e-i and e% are the bulk 
dilational strain, deviatoric strain and shear strain re- 
spectively. The simplification we make is to assume that 



it is adequate to consider the deviatoric strains to be the 
principal order parameter. This is a reasonable assump- 
tion considering the fact that volume changes are usually 
small. Henceforth, we denote the deviatoric strain by e(r) 
and define 

^ = ( d ^- d -^)/V2 = e x (r)-e y ^ (3) 



dx 



dy 



where u x and u y are respectively the displacement fields 
in the x and y directions. 

The free-energy functional of our system with the 
order-parameter e is written as 



F{e(r)} = F L {e(r)} + F LR {e(r)}, 



(4) 



where is a local free-energy functional and Flr is a 
nonlocal long-range term that describes transformation 
induced strain-strain interaction. In a scaled form, we 
write the local free-energy Fl as 



F L = / df 



fl(e(f)) + j(Ve(f)f-a(r)e(r) 



(5) 



where D and a are in a scaled form. The latter represents 
the stress field due to localized defects in the crystal. (In 
a real crystal, lattice defects like dislocations and grain 
boundaries act as sources of stress concentration.) In the 
above equation, /;(e(r)) is the usual Landau polynomial 
for a first order transition given by 

f lW )) = T -e{r) 2 -e{r) A + l -e{rf. (6) 

Here, t — (T ~ T c )/(Tq — T c ) is the scaled temperature. 
To is the first-order transition temperature at which the 
free energy for the product and parent phases are equal, 
and T c is the temperature below which there are only 
two degenerate global minima e = ±cm ■ The stress field 
u{r) in Eqn. 5 modifies the free-energy // in such a way 
that the austenitic phase is locally unstable leading to 
the nucleation of the product phase. 

The physical cause of the long-range interaction is the 
coherency in strain at the parent-product as well as the 
product-product interfaces. An effective long- range in- 
teraction between the deviatoric strains of the trans- 
formed domains has been shown to result from the elim- 
ination of the other strain components, ei and e^, using 
St. Venant compatibility constraint |23j|2J|- Apart from 
determining the dependence on r, one important feature 
is that such a kernel picks out the correct habit plane 
directions which in the present case are [11] and [11]. 
As one of our objectives is to keep the model as simple 
as possible, we have resorted to introducing a long-range 
kernel in a phcnomenological way retaining the feature 
that allows the growth of martensitic domains along the 
habit plane directions 3 5]. (In the present case also, it 
is possible to use this approach to obtain the appropri- 
ate kernel). Wang and Khachaturyan |2^| have shown 
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that the interface can be described by accounting for co- 
herency strains at the parent-product interface by includ- 
ing symmetry allowed fourth order anisotropic long-range 
interaction in the free energy, ie, a term which is invari- 
ant under e — > —e. We define the long-range interaction 
by 

F LR {e} = ~[ J dfd?G{f-P)e 2 {f)e 2 {P). (7) 

As much as in the physical situation, in our model also, 
long- range interaction plays an important role in de- 
scribing the growth of martensite domains. The kernel 



G{f 



is best defined by considering the Fourier rep- 



resentation of the long-range term given by 

F LR {e} = \J dk B (^{e 2 {r)} k {e 2 (r)} k ,, 



(8) 



where {e 2 (f)}k is the Fourier transform of e 2 (r) defined 
as 



{e 2 (r)}k = 



df 



e 2 (f)exp(ik ■ r). 



(9) 



The quantity {e 2 (r)} fc , is the complex conjugate of 
{e 2 (f)} fe . The direction dependent kernel B(k/k) con- 
tains information about the crystallographic details of 
the crystal and defines the habit plane. Apart from the 
favorable directions of growth of the product phase along 
[11] and [11], the free-energy barriers should be large 
along the [10] and [01] directions. These features are 
well captured by the simple kernel 



B 



= -/30(k - A)fc 2 fc 2 



(10) 



where k x and k y are the unit vectors in x and y directions 
(The step function 9(k — A) has been introduced to im- 
pose a cutoff on the range of the long-range interaction.) 
The constant (3 is the strength of the interaction. This 
kernel incorporates the effect of the interface in a simple 
way as the cost of growth progressively increases with 
the transformation in directions where the kernel is pos- 
itive which not only aids growth along the habit plane 
directions but also limits the growth of domains trans- 
verse to it. We stress that this is only a simple choice 
and is not unique. Other kernels with similar orientation 
dependence will give similar results [35j. The real space 
picture of B(k/k) is similar to the long-range interaction 
of Kartha et al [3|. 

Even though we have taken the deviatoric strain as 
the order parameter, basic variables are the displacement 
fields. Thus, we start with the Lagrangian L = T — 
F, where F is the total free-energy and T is the kinetic 
energy associated with the system. The kinetic energy is 
given by 



T = dfp 



du x (f,t) 
dt 



dUy(f, t) 

dt 



(11) 



Here p is the mass density. As mentioned in the introduc- 
tion, since the parent-product interface moves rapidly, 
it is associated with a dissipation which we have repre- 
sented by the Rayleigh dissipative functional |33j . Fur- 
ther, since deviatoric strains are the dominant ones, dis- 
sipative functional is written entirely in terms of e(r), 



R=- 2 i 



d 2 
dr{ Wt e{^t)) . 



(12) 



(Here we have assumed that the bulk and shear strains 
equilibrate rapidly and hence do not contribute to the 
dissipation function.) Now consider the possibility of 
relating the above term to acoustic energy. As far as 
we know there are no attempts to capture the essential 
features of acoustic signals in the context of martensitic 
transformation. To model the AE signals, we recall that 
the mechanism of generation of the AE signals is gener- 
ally attributed to the sudden release of the stored strain 
energy. One area where there has been some efforts to 
model AE signals is in plasticity. In this case, the pro- 
duction of AE signals is attributed to the abrupt motion 
of the dislocations. Consequently, the energy of AE sig- 
nals, E ae (r) is taken to be proportional e 2 (r), where e 
is the local plastic strain rate [3y]. However, in general 
there is spatial inhomogeneity. Then the leading contri- 
bution to total energy E ae cx J(Ve) 2 d 3 r. This clearly 
has the same form as the Rayleigh dissipation functional 
|33j arising from the rapid movement of a localized re- 
gion. Thus, while comparing results of the statistics of 
AE bursts, we need to simply compute R(t). 

We derive the equations of motion for e using the equa- 
tions of motion for the displacement fields given by 

d ( SL\ 5L SR 

i = x,y. (13) 



dt V §Ui 



Su, 



Using the above equation, after computing the functional 
derivatives, we get (see Appendix for details) 



d 2 

'dt 2 



SF 



■7^(r,t) 



(14) 



_8e(f,t) 

which after scaling out p and D can be written in the 
form (in terms of rescaled space and time variables) 



9 2 
dt 



df(r,t) 
de(r, t) 



-cr(f)- V 2 e(r,i) + 7 |-e(r,i) 



at 



+ 2e(r,t) / dkB(k/k){e 2 (k,t)} k e 



ik. 



(15) 



Here, both (3 and 7 are to be taken as rescaled parame- 
ters. The structure of Eqn. (15) is similar to that derived 
in 30] for 1-d except for the long-range term. In the next 
section, we use the above equation to study the morpho- 
logical evolution during martensitic transformations. 



III. NUMERICAL SIMULATIONS 

We now describe the results of our numerical sim- 
ulation for the morphological features. We discretize 



5 



Eqn.(15) on a N x N grid using the Euler's scheme with 
periodic boundary conditions. The mesh size of the grid 
is Ax = 1 and the smallest time step At = 0.002. Most 
results reported here correspond to N — 128, 256. How- 
ever, wherever necessary, we have carried out simulations 
for higher N. A psuedo-spectral technique is employed 
to compute the long-range term |3?| . In this method, we 
compute the discrete Fourier transforms of e 2 (r, t) and 
G(r), take the product and then calculate the inverse 
Fourier transform. In all simulations reported in the pa- 
per, the cutoff A in the long-range expression defined in 
Eq.(10) is chosen to be 0.2. The inhomogeneous stress 
field a(f) is appropriately chosen to describe the defect 
configuration (see below). We consider two situations 
corresponding to the nucleation at a single- defect site 
and at several defect sites. 



A. Nucleation at a Single-Defect Site 

We begin with the study of nucleation and growth of 
domains under a single quench, starting from the austen- 
ite phase to the martensite phase. In general, the stress 
field due to several defects located at fj can be described 
by 

a^=^a Q {r-)exp ^-p ? y (16) 

where ao {fj ) is the magnitude of the stress field at sites 
fj which are randomly chosen defect sites, j max is the 
total number of defect sites, and Q is the width of the 
field. 

We first consider a single isotropic defect with its core 
located at the center of the system. Although, a single 
spherically symmetric defect is a rather artificial system, 
it is useful to clarify the physics of nucleation and growth 
in martensitic systems. For a single-defect case, fj = fo- 
We choose C = 1 and a = 0.3. With this value of ao, the 
system becomes locally unstable at the core r*o. ( The 
threshold value of ao which makes the austenite phase 
unstable at that point is first determined and a slightly 
larger value is used. ) The parameters chosen for the sim- 
ulations are j3 = 50 and 7 = 4. At t — 0, we start with 
e(r, 0) distributed in the interval [—0.005,0.005] repre- 
senting the austenite phase and simultaneously 'turn on' 
the stress field a(f) as we quench the system to r = —2.0. 

Figure 1 shows the nucleation and growth of the 
martensite domains from the defect core at various in- 
stants of time. (Grey regions represent the austenite 
phase with e = 0, black and white regions represent the 
two variants of martensite characterized by e = ±e eg . ) 
We note here that the magnitude of strain in the two 
domains \e e q\ are larger than 6m ( obtained from Eq. 
(6)) due to the contribution arising from the long-range 
term. In a short time after quench (t ~ 7.5), we observe 
the emergence of a small nucleus with e = e eq . In addi- 
tion, we also see the emergence of domains of the other 



variant (e = — e eq ) adjacent to the nucleus in the [11] 
and [ll] directions. (We remark here that the mecha- 
nism of twin ning in this model is the same as that dis- 
cussed in 3Q, :32|, i.e, kinetic energy minimization in the 
presence of dissipation.) The structure further develops 
into twinned arrays, propagating along [11] and [ll] di- 
rections, as can be clearly seen in the snapshot at time 
t = 16.25. This snap shot also reveals an interesting 
feature namely the creation of nuclei ( e — e eq ) located 
close to the +x and — y directions, located at a finite 
distance from the propagating arrays. As the twinned 
fronts propagate, several additional nuclei are created at 
finite distances from the original propagating fronts as 
can seen from the snap shot at t = 18.5. This can be at- 
tributed to the accumulation of long -range stress fields 
at these sites. These new nuclei give birth to secondary 
fronts that also propagate along [11] and [11]. These ob- 
servations in our simulation are in accordance with the 
collective nucleation mechanism discussed in [2^] and the 
experiments by Ferraglio and Mukherjee |14| . The prop- 
agation of these new secondary fronts continues till they 
'collide' with the pre-existing martensite domains and 
stop (t = 20). The morphological evolution eventually 
stops beyond t = 50. 

Figure 2, shows the corresponding evolution of the area 
fraction (f> of the martensitic phase (o). The area fraction 
is computed by counting the number of points on the grid 
for which |e(r, t)\ > 0.5. The transformation is seen to 
start around t ~ 15. The fraction increases sharply till 
about t ~ 20, beyond which it saturates to a value close 
to 0.31. The spurt in the growth between t = 15 and 
20 roughly coincides with the creation of the first set of 
additional nuclei. 



B. Nucleation at Several Defect Sites 

Now we consider a more realistic case of several defect 
sites where the nucleation can occur. Here, we assume a 
random distribution of defects. In the present simulation, 
we choose j m ax = 16 (nearly 0.1% of the total number of 
sites and N =128) and consider <7o(fj) to be uniformly 
distributed in the interval [—0.3,0.3]. All other param- 
eters are same as that for the single-defect case. Ini- 
tially, the system is in a homogeneous state with e(r, 0) 
uniformly distributed in the interval [—0.005,0.005]. At 
t = 0, we turn on the stress-field a(f). In Fig. 3, we 
show the evolution of the system at specifically chosen 
instants of time. As in the case of single- site nucle- 
ation, around t = 15, nucleation of the product phase is 
seen to occur at several sites. By t — 17.5, these nuclei 
grow into twinned lenticular shape. Several additional 
nuclei emerge at finite distance from these original do- 
mains. These new sites at which the product phase nu- 
cleates most often coincides with the pre-existing defect 
sites. However, occasionally, nucleation does occur at 
sites where there were no defects due to stress accumu- 
lation arising from the long-range term, as in the case 
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t= 7.5 t= 16.25 




t = 18.75 t = 50.0 




FIG. 1: Morphological evolution for nucleation at a single- 
defect with /3 — 50, 7 = 4 and r = —2.0. Grey cells corre- 
spond to the austenite phase and the black and the white to 
the two martensite variants. 



t = 15.0 t = 17.5 




FIG. 3: Morphological evolution for multi-defect case with 
defect density 0.1% for /3 = 50, 7 = 4 and r = -2.0 
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FIG. 2: Plot of the transformed area fraction 4>a with respect 
to time for N=128, /3 = 50, 7 = 4, and r = —2.0. o correspond 
to single-site case and • corresponds to the multi-site defect 
case. 

of single-site nucleation. As can be seen from the snap 
shot at t = 15, there is a rapid growth of the product 
phase along [11] and [11] directions forming a cris-cross 
pattern of the martensite domains. We find that there 
is very little growth beyond t — 25 and by t = 50 the 
growth practically stops. A comparison of this final con- 
figuration with the corresponding single-site nucleation 
shows that the final morphologies are very similar im- 
plying that the morphology evolution is independent of 
the original defect ( stress-field) configuration. As for 
the morphology, we can see thin needle-like structures 
emerging from larger domains. Figure 2 also shows the 
time evolution of the area fraction cf> for multi-defect case 
(•) which is very similar to the single-defect case, which 



again emphasizes the fact that the time evolution is not 
sensitive to initial defect configuration. All the adaptive 
domain adjustments take place beyond this time regime. 
Thereafter, the area fraction saturates. 



We briefly remark on the influence of the two param- 
eters (3 and 7 on the morphology. The quantity 7 repre- 
sents the strength of dissipation and (3 the strength of the 
transformation induced long-range interaction. We find 
that the lateral width of the arms for smaller (3 is larger 
than that corresponding to larger values of (3 for a fixed 
value of 7. The area fraction is higher for smaller [3 which 
can be attributed to the fact that a lower value of (3 cor- 
responds to a lower interaction between the transformed 
regions and hence, lower energy cost to grow transverse 
to the direction of propagation. ( Note that (3 = corre- 
sponds to symmetric growth.) Now consider the influence 
of the dissipative term on the morphology. Our finding is 
that the twin width is larger for larger values of 7. This 
result can be understood by noting that the overdampcd 
case (7 >> 1) corresponds to growth as a single variant 
[30|. With the inertial effects, growth as a single variant 
is prohibited by high kinetic energy cost. In view of this, 
for large values of 7, the dissipative term dominates and 
the system grows as a single variant for a larger distance 
than that for the low damping case. 
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FIG. 4: Area fraction of the transformed phase <j>A as a func- 
tion of r for both single defect (o) and multi-defect (•) nu- 
cleation cases. The parameter values are 13 = 50, 7 = 4 and 
N = 128. The defect density in multi-defect is 0.1% of sites. 



IV. THERMAL CYCLING AND HYSTERESIS 
A. Morphological Features 

One key feature of a martensitic transformation is the 
hysteresis observed when the system is subjected to ther- 
mal cycling. Here, the system is cooled at a specific rate 
from the high temperature phase to the low temperature 
phase and then subsequently heated back to the high 
temperature phase. We have performed 'continuous cool- 
ing' and 'heating' computer simulations where we change 
t at a constant rate: the interval t = 40 to -80 is cooled 
in 1000 time steps. We have monitored both the mor- 
phology and the area fraction of the transformed phase 
for single-defect nucleation as well as multi-defect nucle- 
ation case with several system sizes ( N=128,256 ). Wc 
use the same initial conditions for both these cases as that 
used for a single quench situation. The initial condition 
for the reverse transformation is the final configuration 
obtained during the cooling run. 

In Fig. 4, we have shown the variation of the area 
fraction <\>a with r for the heating and the cooling runs 
for the single defect case (•) and multi defect case (o) 
for N = 128. In the cooling run, for the multi-defect 
case, the transformation starts around r ~ —2.0 showing 
a rapid increase in <j>A ( ~ 30%). Thereafter, there is 
a nearly linear increase up to 90% at which the growth 
rate tapers off. The system gets fully transformed at 
T m f ~ —60. In the heating run, the reverse transforma- 
tion does not start till r as <~ —22.0 and thereafter, <f>A 
decreases almost linearly till the transformation is nearly 
complete around r <~ 18. As can be seen from the fig- 
ure (Fig. 4), the difference between the hysteresis cycles 
corresponding to the multi- site defect nucleation and 
single-defect nucleation cases is small. 

We have followed the morphological evolution of the 
martensitic domains both for the single and multi-defect 
cases. Here, we shall only discuss the more realistic multi- 
defect case. Figure 5 shows the snap shots of the pat- 
tern at specifically chosen values of r as the system is 
taken through a cycle. In the snapshot corresponding 



x = -3.5 x = -6.5 




x = - 24.5 x = - 45.5 




FIG. 5: Morphological evolution during a cooling and heating 
cycle. The parameter values are (5 = 50, 7 = 4 and N = 128 
with a defect density 0.1% of the sites. 



to t = —3.5, one can see the nucleation of martensitic 
domains at multiple locations. As the system is further 
'undercooled' to lower values of t, not only does these 
twins propagate in the [11] and [11] directions, the thick- 
ness of the martensite domains also increases. This can 
be seen from the snapshots corresponding to r = —6.5 
and t = — 24.5. As the system is further undercooled, 
the twin width further increases with the last regions 
transforming around r = —59. A late stage snap shot at 
t = —45.5 corresponding to cf) <~ 90% is shown. For the 
heating cycle, the final configuration of the cooling run 
is taken as the initial configuration. The austenite phase 
appears around r <~ —22. As r is further increased, the 
martensite phase can be seen to be gradually disappear 
in the snapshots corresponding to r = —5 and 13. Fur- 
ther, one can see that the overall morphology in the final 
stages of the heating run is significantly different from 
the initial stages of the cooling run. Thus, in our model, 
there is no long term memory effects, though there is 
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short term memory as will be shown later. Finally, by 
t = 18.0, the martensite phase disappears completely. 



B. Power Law Statistics During Thermal Cycling 

From Fig. 4, it appears that changes in (pA are smooth 
on the scale shown in the figure. In reality, on a finer 
scale, the changes in (f)A are actually jerky. In fact, in ex- 
periments, thermal cycling is accompanied by the emis- 
sion of acoustic energy in the form of bursts, a feature 
that reflects the jerky nature of the transformation. In 
the model, as mentioned earlier, the energy of acoustic 
signals is captured by the rate of energy dissipated given 
by R(t) = —dE/dt. We have calculated R(t) during the 
heating and cooling runs. Figure 6 shows R(t) as a func- 
tion of r with the inset showing the enlarged section of 
the peak. The figure clearly shows that the rate of energy 
release occurs in bursts consistent with acoustic emission 
studies [j| during thermal cycling. 

Since, in experiments one finds that the AE signals 
show a power law statistics, we have investigated the dis- 
tributions of the amplitudes of the AE signals and their 
durations. Denoting the amplitude of R(t) by Ra, we 
find that the distribution D(Ra) of Ra has a tendency 
to approach a power law, ie., D(Ra) ~ A~ aR with an ex- 
ponent otR. Figure shows a log-log plot of D(Ra) as a 
function of Ra, for both the single-site nucleation case (•) 
and the multi-defect case (o). From the figure, it is clear 
that both these cases exhibit the same exponent value 
an ~ 2.5 over three orders in D(Ra)- In experiments, 
one also finds that duration of these bursts also obey a 
power law statistics. To verify this, we have also plot- 
ted the distribution D(At) of the durations At of energy 
bursts for both the single and multi defect cases. We find 
that D(At) ~ At~ TR with an exponent value tr w 3.2, 
although, the scaling regime is not as impressive as for 
Ra- ( Here, we remark that typically the scaling regime 
for the durations of the events is much smaller than that 
for the amplitudes even in models of SOC H, l3S|.) We 
have also calculated the conditional average < Ra > c 
for a given value of At |31|. This is expected to obey a 
power law given by < Ra > c ~ At XR . The value we get 
is about x w 1.36 for both these single and multi defect 
nucleation cases. Using these values, we find that the 
scaling relation a — x(j3 — 1) + 1 is satisfied quite well. 
We have also carried out a similar analysis on R(t) for the 
heating run. Even though, the changes in R(t) occurs in 
bursts, we find the scatter is considerably more than for 
the cooling run. Wc will comment on this aspect later. 

A comment is in order regarding the exponents. In 
experiments one actually plots the distribution for the 
amplitude of the AE signals, while in our model acous- 
tic energy plays a natural role. Thus, while compar- 
ing the exponents values, one needs to use Ra ~ A AE , 
where Aae is the amplitude of the AE signal. Using 
the relation between the two joint probability distribu- 
tions D(RA,At) oc D(Aae, At)/AAE: one easily finds 
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FIG. 6: R(t) as a function of r with inset showing enlarged 
section of the peak during cooling. The parameter values are 
13 = 50,7 = 4 and N = 256. The defect density in multi- 
defect is 1% of sites. 
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FIG. 7: log-log plot of D(R A ) as a function of R A . The 
parameter values are /3 = 50, 7 = 4 and N = 256. The defect 
density in multi-defect is 1% of sites, [o corresponds to single- 
defect case and • corresponds to multi-defect case.] 



that an = [olae + l)/2 with the other two exponent val- 
ues remaining unchanged. Using the experimental val- 
ues of olae ~ 3.8, (jae ~ 3.8 and xae ~ 1), we see 
that the calculated value of an « 2.4. Considering the 
fact that our model is two dimensional, we see that the 
agreement of the exponent values is reasonable. Here, it 
should be pointed out that even the number of marten- 
site variants in 3-D are more than in 2-D. Thus, it would 
be unrealistic to expect that the mechanisms in 3-D can 
be accounted for in 2-D. 

We now explain the origin of the power law statistics in 
the model. This can be traced to the fact that we have in- 
cluded important ingredients of SOC dynamics, namely, 
the threshold dynamics, dissipation, the generation of 
large number of metastable states during the transfor- 
mation, and a relaxation mechanism for the stored en- 
ergy. The relaxation of the stored energy occurs at very 
fast time scales comparable to time scale of the speed of 
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et al have verified that in martensite transformation, it is 
the dynamical disorder ( transformation induced) that is 
at the root of the avalanches. This exactly what is well 
captured by the present model. 



SHAPE MEMORY EFFECT 



A. Correlated behavior of AE signals 



FIG. 8: (a). Area fraction cj> during thermal cycling in the 
interval r = -8 to 4 for cycles 1 (o),2 (A), 3 (»),4 (>). 
(b)Hysteresis for the full cycle from the austenite to the 
martensite phase and back. (• ) is the starting point of the 
small thermal cycles. 



sound as can be seen from the fact that the basic vari- 
ables are the displacement fields. ( We note that this is 
the fastest time scale.) Indeed, from our simulations we 
find that the interface movement occurs at time scales of 
a few units of (scaled) time as can be inferred by rapid 
increase in the area fraction of the transformed phase, 
typically of the order of a few units of scaled time. (See 
Fig. |2) Compared to this the driving force generated by 
thermal cycling increases with temperature slowly which 
is one of the characteristic features of SOC dynamics. 

Another important feature of the model, as also that 
of SOC dynamics, is the creation of large number 
metastable states during cooling or heating runs. This 
is a direct consequence of an interplay between the local 
free energy ( free energy barrier) and the long-range inter- 
action between the transformed domains as can be seen 
from the following reasoning. We note that the value of 
the long-range term at any spatial location is the result 
of the superposition of the contributions arising from the 
spatial distribution of the already transformed domains. 
As a consequence, the free energy surface at any given 
time is a complex terrain of local barriers ( metastable 
states). It must be noted that these local thresholds are 
self generated ( transformation induced). At a given time, 
these local thresholds must be overcome by the increase 
in the driving force arising from the slow cooling (or 
heating). We note that once a local barrier is overcome, 
part of the driving force goes in creating a new twin and 
the rest is dissipated in the form of burst of energy due 
to the advancing one or more interfaces. The fact that 
long-range interaction is at the root of creating the local 
thresholds is further supported by the fact that we find a 
power law distributions even in the single- site nucleation 
case. (See • in Fig. 0) The presence of defect sites only 
serves to trigger the initial nucleation process. This must 
be contrasted with disorder based Ising models 25] which 
also produce power law statistics for avalanches and field 
induced hysteresis. However, as mentioned earlier, Vives 



As mentioned in the introduction, one experimental 
result (known for some time) is the highly correlated be- 
havior of AE signals when the system is subjected to 
thermal cycling in a small temperature interval. To ver- 
ify if this result can be captured by our model, consider 
the system being cooled from the austenite phase to a 
point in the martensite phase where a desired amount of 
martensite phase has developed. Starting from an appro- 
priate point ( shown as • in Fig. |8Jd, here chosen to be 
4> ~ 0.58 in the full thermal cycle), we subject the sys- 
tem to repeated thermal cycling in a small temperature 
range r m i„ = —8 and r max = 4. As in the case of full 
thermal cycling, for the small heating cycle also, the final 
configuration attained at r max — 4, is taken as the initial 
configuration for the cooling run. Calculations have been 
performed for a range of parameter values of 50 < (3 < 10 
and 5 < 7 < 1. For the present calculation, we have used 
f3 = 35 and 7 = 4. The value of used here is higher 
than that reported in our earlier paper . We shall dis- 
cuss the influence of increasing (3 soon. Other parameter 
values are A = 0.2, Q = 1 and N = 128. During the first 
few cycles, the loops in the area fraction <j> verses r drift 
slightly, but stabilize after a first few cycles, here, after 
the sixth cycle. The first few cycles play the role of the 
training period known in experiments. After these first 
few cycles the system eventually circulates in the same 
set of configurations as we will show. ( The number of 
training cycles in the model is only few, however in ex- 
periments, the training period is typically a few cycles 
for CuAlZn ,9j, but could be much larger in some other 
alloys.) 

During the training cycles, the energy dissipated R(t) 
evolves continuously, stabilizing only after the training 
period. A plot of R{t) for several forward and reverse 
cycles ( seventh to tenth ) after stabilization is shown in 
Fig. It is clear that the energy bursts (which mimic 
the AE signals), as in experiments |39l l40l] . exhibit a near 
repetitive pattern in time (temperature) during succes- 
sive heating and cooling parts of the cycles. As can be 
seen the bursts are much more noisy when compared to 
slightly smaller value of (3 used in our earlier study [T]|. 
Further, we find that increasing f3 requires more cycles 
to stabilize ( for instance, the cycles stabilize after the 
fourth for f3 = 25). 
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FIG. 9: Repetitive nature of R(t) for cycles 7 to 10. 



B. Connection to shape memory effect 

In order to establish a correspondence between R(t) 
and the changes in the spatial configuration of the 
martensite domains, we have simultaneously monitored 
the morphology over all the cycles. We find that the 
morphology changes drastically during the first few cy- 
cles even though the macroscopic state of system in terms 
of <j) — t returns to nearly the same point at the end of 
each cycle. The domain configurations at the beginning 
of first few cycles are shown in Fig. Figure EH a is 
the starting morphology for the first cycle that has been 
obtained by slowly cooling from the austenite phase. Fig- 
ure HOB is that obtained at the end of sixth cycle. As can 
be seen, most changes occur during the first few cycles 
with large number of changes occurring in the first cycle 
itself. We find that during the first cycle itself, most of 
the curved twin interfaces in the initial configuration ( 
Fig. 1 10b ) are rendered straight and several small twins 
coalesce to form a single variant of the same type. One 
can also notice that some regions of the austenite phase 
separating the martensite domains are washed out. Sub- 
sequent cycles also have the same effect, but less effective. 
The morphology stabilizes even though some curved in- 
terfaces with small regions of martensite phases remain ( 
see the top left edge in Fig. HUH ) . After the sixth cycles 
very little changes could be detected. 

Now consider the changes in the morphology during a 
single cycle in the region where R(t) is repetitive. The 
snapshots of the morphology during one such stabilized 
cooling and heating cycle, the seventh one, (starting from 
the initial configuration shown in Fig. I1UH ) at selected 
intervals is displayed in Fig. ^] It is clear that during 
heating, the martensite domains shrink, opening up the 
austenite phase and some martensite domains even dis- 
appear. However, during cooling these domains reappear 
and the eventual morphology at the end of the cycle is 
practically recovered on returning to the starting point 
on the (<fi, t) diagram. As can be seen, the final configu- 
ration obtained during the seventh cycle, Fig. can be 
seen to be practically the same as Fig. \1(A d which is the 
initial configuration for the seventh cycle. These obser- 
vations are consistent with that observed in experiments 
□ 




FIG. 10: Morphology of the initial configurations at r = —8.0 
for the first, second, fourth and the sixth thermal cycles. 



Now, we attempt to provide a physical explanation for 
the repetitive nature of the energy bursts and the reversal 
of the morphology under thermal cycling. To understand 
this, we need to understand the role played by the train- 
ing cycles. Indeed, the repetitive pattern of the energy 
bursts during successive cycles (after the training period) 
is an indication that the system traverses through the 
same set of metastable states. From Fig. 110b ,. we note 
that the initial configuration for the first cycle has a large 
number of small domains compared to the configuration 
obtained after the sixth cycle shown in Fig. ^|d. In ad- 
dition, the twin interfaces of Fig. HUfa are rough. ( The 
interface is considerably more curved for small [3 values. 
Compare Fig. 3 of Ref. [TT|| . For the present case, it is 
relatively straight.) Such configurations are generally ex- 
pected to have higher energy compared to straighter ones. 
Thus, the initial configuration used for cycling (Fig. II Oh ) 
corresponds only to a local shallow minimum. During the 
first few cycles, the free energy landscape is so modified 
that it smoothen's out the energy barriers corresponding 
to large number of twin interfaces in Fig. HOfa ) with very 
little change in the area fraction. As in the case of power 
law distribution, a crucial role in smoothening process is 
actually played by the long-range interaction term, as the 
growth (shrinkage) of a martensite domain is influenced 
by the configuration of rest of the domains. To verify 
this, we have computed the free energy Flr arising from 
the long-range interaction between the domains and find 
that it actually becomes more negative with successive 
cycles saturating after first few cycles. This additional 
contribution leads to a reduction in the local free energy, 
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FIG. 11: Sequential morphological snapshots for r = 
1.6(o), 4.0(&), -2.0(c), -8.0(d) during the seventh cycle. The 
initial configuration for the cycle is that shown in Fig. 10(d). 

Fl, as well. The net effect is to create a deeper set of 
metastable states for the system to circulate for the sta- 
bilized cycles. Within one such stabilized cycles, say sev- 
enth, the starting configuration (Fig. llOb O has the lowest 
free energy reaching a maximum at the end of a heating 
cycle, ie., at r = 4.0, Fig. II lb . Thereafter it decreases 
during cooling. Note that the morphology at t = 4.0 is 
significantly different from the starting morphology. The 
increasing or decreasing width of the martensite domains 
as we decrease or increase the temperature is surprisingly 
similar to that observed in experiments. ( See Fig. 12 of 
Ref.0.) 

VI. PRETRANSITIONAL EFFECTS 

Now we consider the possibility of recovering pretran- 
sitional effects in our model. As mentioned in the intro- 
duction, the mechanism attributed by Kartha et al is the 
dependence of the elastic constant on local disorder in 
composition. It is clear that once the basic mechanism 
is included, the model should lead to results similar to 
that in Kartha et al. We shall adopt the same idea and 
assume that the transition temperature depends on local 
compositional fluctuations through 

T (f)=T (c)-ASc(f), (17) 

where A is the relative strength of coupling to composi- 
tional fluctuations 5c, assumed here to be randomly dis- 
tributed. In scaled variable that we use, this leads to r 




FIG. 12: Development of Tweed like pattern 



being replaced by r = r(c) + adc, where c is the average 
concentration and a is a scaled variable. 

In our simulations, we assume that the initial concen- 
tration fluctuations, considered as frozen, are drawn from 
a normalized Gaussian. In our case, we just need to solve 
the equations of motion numerically for various values of 
r and a. (We have varied a from 0.0 to 1.2.) The range 
of values of (3 and 7 wherein we find the tweed struc- 
ture is in the region of relatively low values. ( (3 < 20 and 
7 < 0.7.) Here, we report results for a — 0.9,(3 = 10 
and 7 = 0.1. A typical set of morphologies are shown 
in Fig. E| The pretransitional effects become notice- 
able even when r = 4. For r = 2 the directionality of 
the tweed pattern is already evident. As we decrease 
the temperature, this structure becomes more dominant. 
It must be mentioned here that in our case, these pat- 
terns do not change after reaching a steady state. Thus, 
we do not have the dynamic tweed structure reported 
by Kartha et al. This is because, in our model the ef- 
fect of temperature goes only in the local free energy 
and there are no thermal fluctuations. We also note that 
the pretransitional effects are pronounced when both (3 
and 7 are small. We shall comment on this later. How- 
ever, we note that including the changes in the transition 
temperature on local disorder is somewhat similar to the 
contribution from defects except that internal stresses ( 
arising from deviations from the average concentration) 
are of much smaller magnitude compared to that due to 
defects. Thus, it is clear that one can mimic the present 
situation through Eq. Ilrjl bv appropriate choice of 00 and 
£j. In our case, the a — r plane, the boundary between 
the austenite and the tweed phases is linear in the range 
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(0,0) to (1.2,4). 



VII. SUMMARY AND DISCUSSIONS 

In summary, we have presented the results of a compre- 
hensive study of the dynamics of strain driven marten- 
sitic transformations within the framework of a two di- 
mensional square-to-rectangle transition. Due to the fact 
that the long-range term is introduced phenomenologi- 
cally, the model should be viewed as a model whose pri- 
mary aim is to capture the essential physics of the trans- 
formations. We however note that it includes all the im- 
portant contributions arising from different mechanisms 
in a transparent way. In doing so, we are able to study the 
dynamics of the strain driven transition that explains the 
three unusual features of the martensitic transformation 
that were sought to modelled. The hrst important feature 
that emerges from the model is the fact that the elastic 
energy stored is released in the form of bursts. (Note 
that we have also established a correspondence between 
the dissipative functional and the energy of AE signals.) 

Second, the model provides a proper basis to explain 
the power law statistics of the AE signals observed in 
experiments. As stated earlier, the power law statistics 
arises due to the fact that we have included threshold 
dynamics, dissipation, the generation of large number of 
metastable states during the transformation, and a relax- 
ation mechanism for the stored energy whose time scale 
is much faster than the time scale of the drive force. It is 
interesting to point out that both cases of single-defect 
and multi-defect nucleation lead to power laws strongly 
suggesting that quenched defects do not play any role in 
the power law. The model also predicts the near repet- 
itive bursts of energy under successive thermal cycles in 
a small temperature interval as observed in experiments 
on acoustic emission 0, |3|| H3 . This comes as a surprise 
as the nature of correlation in the latter case is 'periodic' 
in contrast to the scale free nature of correlations dis- 
cussed in the power law statistics of AE signals. More 
importantly, these bursts of energy have been shown to 
be correlated to the growth and shrinkage of martensite 
plates. The underlying cause of the correlated nature of 
AE signals and reversal of the morphology during suc- 
cessive cycles after the training cycles has been traced to 
the influence of the long-range interaction term. The role 
of training cycles is also elucidated. During the training 
period the long-range term has a tendency to smoothen 
out higher energy barriers in the free energy landscape. 
This in turn induces a transformation pathway along a 
unique set of low energy metastable configurations. We 
expect that our analysis provides a good insight into the 
shape memory effect which finds immense applications 
in a variety of areas from mechanical actuators to bio- 
medical applications ^2 • To the best of authors knowl- 
edge, this is the first model which shows near full reversal 
of morphology under thermal cycling. Finally, it is not 
surprising that the model also reproduces tweed like fea- 



tures. 

A comment may be in order on the physical interpre- 
tation of the cause of the memory effect in our model. 
Recall that during the course of the transformation (ei- 
ther in single quench situation or in thermal cycling), 
nucleation occurs at sites that are not necessarily the 
sites of quenched defects, but arising from constructive 
superposition of long- range stress fields of the preexist- 
ing domains. However, as quenched defects have been 
represented by a Gaussian stress profile in our model, in 
the same spirit, these new nucleation sites could in prin- 
ciple be interpreted as creating new defects at these sites. 
(See the comment in Ref. |41|.) We also believe that it 
may be possible to model this with an additional degree 
of freedom for the defect kinetics. We expect that our 
analysis provides a good insight into the shape memory 
effect which finds immense applications in a variety of 
areas from mechanical actuators to bio-medical applica- 
tions 0. 

It is interesting to note the similarity of morphological 
patterns with real micrographs. Patterns studied in the 
context of single-defect and multi-defect nucleation show 
that the eventual morphology are independent of the the 
original defect configuration. It is not clear if this is 
true of real systems although experiments of Vives et al 
provide an indirect evidence. In addition, the model 
shows hysteresis under thermal cycling. Even though, 
the model explains several generic features stated above, 
we stress that our model is not material specific. 

As mentioned the pretransitional effect is pronounced 
when both (3 and 7 values are small. The low values of 
7 is physically understandable as one expects that there 
would be hardly any dissipation in the high tempera- 
ture phase where there are no martensite domains. The 
smallness of (3 is also understandable physically (as the 
magnitude of fluctuations in the composition is not high). 
However, the value of f3 cannot be determined in terms 
of the elastic constants as in the case where the kernel is 
derived by using the compatibility relation Q . 

Some comments may be in order on SOC type of fea- 
tures obtained in the model. We note here that unlike 
most SOC models where noise is essential 0, |38|, the 
model is fully dynamical in the sense that noise has no 
role in the generation of power law statistics. In spirit, 
the model is closer to that by Gil and Sornette^l where 
an explicit threshold term is introduced in the form of a 
subcritical Hopf bifurcation. However, it must be pointed 
out that noise is essential in their model as well. Lastly, 
to the best of the authors knowledge, this is one of the 
few fully continuous space-time dynamical (noise free) 
model for SOC. 
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APPENDIX 



and a similar equation for u y . One also finds 



Here we present some details of the derivation of Eqn. 
15. To start with we evaluate , 
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Using the definition for e, the above expression reduces 
to a compact form 
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An analogous calculation for holds except for a neg- 
ative sign on the left hand side with interchange of the x 
with y. Similarly, one can show that 
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and a similar equation for u y . The functional derivative 
(x component) for the dissipation functional is 
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Then, using these results in Eq.(13), and taking appropri- 
ate derivatives with respect x and y, we get the equation 
for the strain order parameter as 



p—*e(r,t) = V 2 
dt 2 v ' 
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Using F — Fl + Fnl, in a compact form the above equa- 
tion can be written as 
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